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1. Introduction 

If the standard paradigm of cosmology is correct, the observed structure of the universe origi¬ 
nated from tiny density fluctuations via gravitational instability. The clustering process is inevitably 
associated with peculiar motions of matter, i.e. deviations from a pure Hubble flow, which exhibit 
a coherent pattern on large scales. Since galaxies can, to good approximation, be treated as test 
particles, they should appropriately reflect the underlying peculiar velocity field which contains 
valuable information for constraining and discriminating between different cosmological models. 

Redshifts of galaxies are systematically altered by the line-of-sight components of their pe¬ 
culiar velocities and differ from their actual distances. Consequently, intrinsic galaxy luminosities 
inferred from the observed flux using redshifts rather than distances appear brighter or dimmer. 
Since this effect is obscured by the natural scatter in the distribution of magnitudes, it cannot be 
used to derive the peculiar velocities of individual galaxies. However, it is possible to approach 
peculiar velocities in a statistical sense by constraining the parameters of some appropriate, prede¬ 
fined velocity model. For instance, constraints on the bulk motion of galaxies in a subvolume of a 
given survey can be derived by comparing the luminosity distribution of galaxies in the subvolume 
with that of the whole survey. Reconstructing the linear velocity field from the observed density 
field in redshift space, this technique further provides a way of estimating the growth rate of density 
perturbations which is independent from the apparent clustering anisotropy of galaxies [|I], [|]. 

Below I briefly review the luminosity-based approach outlined above and discuss its recent 
application to galaxy data from the Sloan Digital Sky Survey (SDSS) [[|]. Methods of this kind 
have first been adopted to estimate the motion of Virgo relative to the Local Group [[jj], and more 
recently, to constrain bulk flows and the growth rate in the local universe within z ~ 0.01 [^|. [hj]. 


2. Velocity-induced modulation of observed galaxy luminosities 


Due to inhomogeneities, the observed redshift z of a galaxy is generally different from its 
cosmological one, z c , which is defined for the homogeneous background. In linear perturbation 
theory, the two quantities are connected through [jT] 

Z-Zc = y(t,r) _ <1>(?,r) _ 2_ r ‘0 d®[rr(t),t] _ V(t,r) 

1 +Z c C 2 C 2 Jt[r) dt C 

where r denotes a unit vector along the line of sight to the galaxy, V is the galaxy’s (physical) radial 
peculiar velocity, and the usual gravitational potential. Focusing on low redshifts at the present 
time, the velocity V is explicitly assumed as the dominant contribution. All fields are considered 
relative to their present-day values at to- The difference in redshifts enters the distance modulus 
DM = 25 + 51og 10 [D/,/Mpc] and causes the observed absolute magnitudes M to deviate from their 
true values M- l \ i.e. 


M = m- DM(z) - K{z) + Q(z) = M (t ) + 5 log 10 


D L (Zc) 
Dl{z) ’ 


where m is the apparent magnitude, Di is the luminosity distance, and the functions Q(z) and K(z) 
account for luminosity evolution and A'-correction [Q]. respectively. Restricted to scales where 
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linear theory is valid, M — M'f' varies systematically over the sky, and thus provides a probe of the 
cosmic velocity field. 

Given a galaxy survey with magnitudes, (spectroscopic) redshifts, and angular positions r, on 
the sky, the idea is to choose an appropriate parameterized model of V (r, z) and to maximize the 
probability of the data, 


Ptot = n PiMilzuViruZi)) = n ( / J2 HM)dM j , 


where redshift errors are neglected [||]. Here <j)(M) denotes the galaxy luminosity function (LF), 
and the limiting magnitudes M ± depend on the cosmological redshift z c , and hence on V ( r,z ). The 
motivation of this approach is to find those velocity model parameters which minimize the spread 
in the observed magnitudes. 


3. Constraints from SDSS galaxy luminosities at z ~ 0.1 

Galaxy data from the SDSS Data Release 7 [S] trace the cosmic velocity field out to z ~ 0.1. 
Here I summarize recent results obtained from applying the luminosity method to suitable subsets 
comprising up to half a million galaxies [ |i~C| . |TT| |. 

Data . —The analysis is based on the latest version of the NYU Value-Added Galaxy Cat¬ 
alog (NYU-VAGC), adopting the subsample safe to minimize incompleteness and systemat- 
ics & Using Petrosian 01 r-band magnitudes, only galaxies with 14.5 < m r < 17.6, —22.5 < 
M r — 5 log 1 0 /r < —17.0, and 0.02 < z < 0.22 (relative to the CMB frame) are included. In addition, 
galaxy mock catalogs were built to study errors and known systematics of the data. 

Radial velocity model. —Constraints on velocity moments and derived quantities assume a bin- 
averaged model V (r) in two redshift bins, 0.02 < z < 0.07 and 0.07 < z < 0.22. For each bin, the 
velocity field was decomposed into spherical harmonics, i.e. 

aim = [ dnv (r)Y lm if) , V(r) = £ a lm Y,* m (r), l > 0, 

^ l,m 

where the sum over l is cut at some maximum value / max . Because the SDSS data cover only 
part of the sky, the inferred a/ m are not statistically independent. The impact of the angular mask 
was studied with the help of suitable galaxy mock catalogs. The monopole term (/ = 0) was not 
included since it is degenerate with an overall shift of magnitudes. 

LF estimators . —Reliable measurements of the galaxy LF form a key step in the analysis. To 
assess the robustness of results with respect to different LF models, the data were examined using 
LF estimators based on a Schechter form and a more flexible spline-based model, together with 
several combinations and variations thereof. For simplicity, a linear dependence of the luminosity 
evolution with redshift was assumed. 

Bulk flows and higher-order velocity moments. —Accounting for known systematic errors in 
the SDSS photometry, measurements of “bulk flows” are consistent with a standard ACDM cos¬ 
mology at a l-2a confidence level in both redshift bins. A joint analysis of the corresponding 
three Cartesian components confirms this result. To characterize higher-order moments, direct 
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Figure 1: Estimates of Cg from galaxy mocks: shown are the recovered distributions and respective Gaussian 
fits with (solid lines) and without (dashed lines) the inclusion of a (randomly oriented) photometric tilt, using 
the information in both redshift bins (left) and the low-z bin only (right). 


constraints on the angular velocity power spectrum C/ = (|a/ m | 2 ) were obtained up to the octupole 
contribution. The estimated C/ were found compatible with the theoretical power spectra of the 
ACDM cosmology. 

Constraints on a*.—Assuming a prior on the C/ as dictated by the ACDM model with fixed 
Hubble constant and density parameters, the amplitude of the velocity field was estimated in terms 
of o%. Because of a dipole-like tilt in the galaxy magnitudes |[T3||, raw estimates of a x were biased 
toward larger values (see Figure |T|). After correcting for this magnitude tilt, it was found that 
Og = 1.1 ± 0.4 for the combination of both redshift bins and as = 1.0 ± 0.5 for the low-z bin only, 
where the low accuracy is due to the limited number of galaxies. 

The linear growth rate. —Modeling the linear velocity field from the observed galaxy cluster¬ 
ing in redshift space, the luminosity approach is capable of constraining the growth rate of density 
perturbations, /3 = f(Q.)/b, where b is the linear galaxy bias [||]. To this end, both magnitude- and 
volume-limited subsamples were selected over a rectangular patch (in survey coordinates) within 
the range 0.06 < z < 0.12. Following [[i~4|l, the velocity reconstruction was done in “harmonic” 
space by smoothing the density field with a Gaussian kernel of 10/i 1 Mpc radius and solving 

J__d / 2 d<F/,„ \ _ 1 /(/+ !)<£>/,„ _ /3 / g _ dlog S d<F/ m \ 

s 2 ds\ ds / 1 + [3 s 2 1 + /3 \ lm ds d s ) ' 


where S is the selection function and <F(.v) is the velocity potential expressed in redshift space. 
Boundary conditions were fixed by setting the density contrast outside the observed volume to 
zero. An example of how the full velocity field modulates galaxy magnitudes at z = 0.1 is depicted 
in Figure || (left panel). The velocity reconstructions and the method’s accuracy were assessed 
with the help of mocks generated from the Millennium Simulation [15, 16]. Excluding the dipole 
in the velocity reconstruction (Z > 1), the found distribution of /3-estimates peaks at the true value 
/3 trU e = 0.52, deviating from a symmetric Gaussian mainly because of the velocity field’s nonlinear 
dependence on /J (see right panel of Figure Q). A preliminary analysis of the SDSS data with fixed 
Hubble constant and density parameters from Jl7[ ] yields /3 = 0.42 ± 0.14 which, using the power 
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Figure 2: Left panel. —Velocity-induced modulation of galaxy magnitudes at z = 0.1 over an angular patch 
expressed in SDSS survey coordinates (based on a random mock). Right panel. —Distribution of /j- estimates 
(and Gaussian fit) derived from mock catalogs, excluding the dipole in the velocity reconstruction {l > 1). 


spectrum amplitude of L'-galaxies from [18], translates into /as = 0.37 ±0.13 and is consistent 
with the ACDM model. For samples with similar characteristics, the accuracy is comparable to 
what is obtained from recent measurements of anisotropic clustering in redshift space [|I9|]. 


4. Outlook 


Current and next-generation spectroscopic galaxy surveys are designed to reduce data-inherent 
systematics because of larger sky coverage and improved photometric calibration in ground- and 
space-based experiments [|20|,pl]]. Together with the above results, these observational perspectives 
give confidence that the luminosity-based approach will be established as a standard cosmological 
probe, independent from and complementary to the more traditional ones based on galaxy cluster¬ 
ing, gravitational lensing and redshift-space distortions. The method considered here does not re¬ 
quire accurate redshifts and can also be used with photometric redshift surveys such as the 2MASS 


Photometric Redshift catalog (2MPZ) | [22| | to recover signals on scales larger than the spread of the 
redshift error. 


References 

[1] N. Kaiser, MNRAS, 227, 1 (1987). 

[2] A.J.S. Hamilton, Linear Redshift Distortions: a Review, in The Evolving Universe, Astrophysics and 
Space Science Library, 231 (1998) [astro-ph/ 97 08102]. 

[3] D.G. York, J. Adelman, J.E. Anderson, Jr., et al.,AJ, 120, 1579 (2000) [astro-ph/000 6396]. 

[4] G.A. Tammann, A. Yahil, & A. Sandage, ApJ, 234, 775 (1979). 

[5] A. Nusser, E. Branchini, & M. Davis, ApJ, 735, 77 (2011) [arXiv: 1102.418 9]. 

[6] A. Nusser, E. Branchini, & M. Davis, ApJ, 744, 193 (2012) [arXiv: 110 6.6145]. 

[7] R.K. Sachs, & A.M. Wolfe, ApJ, 147, 73 (1967). 


5 




















The linear velocity field ofSDSS DR7 galaxies 


Martin Feix 


[8] M.R. Blanton, & S. Roweis, AJ, 133, 734 (2007) [astro-ph/0 60 6170]. 

[9] K.N. Abazajian, J.K. Adelman-McCarthy, M.A. Agueros, M.A., et al.,ApJS, 182 543 (2009) 

[arXiv: 0 812.0 64 9]. 

[10] M. Feix, A. Nusser, & E. Branchini, JCAP, 09, 019 (2014) [arXiv: 1405.6710]. 

[11] M. Feix, A. Nusser, & E. Branchini, in preparation. 

[12] M.R. Blanton, D.J. Schlegel, M.A. Strauss, et al.,AJ, 129, 2562 (2005) [astro-ph/0410166]. 

[13] N. Padmanabhan, D.J. Schlegel, D.P. Finkbeiner, et al.,ApJ , 674, 1217 (2008) 

[astro-ph/07 03454], 

[14] A. Nusser, & M. Davis, ApJ (Letters), 421, LI (1994) [astro-ph/9309009]. 

[15] V. Springel, S.D.M. White, A. Jenkins, et al., Nature , 435, 7042 (2005) [astro-ph / 05 04 0 97]. 

[16] B.M.B. Henriques, S.D.M. White, G. Lemson, et al., MNRAS. 421, 4 (2012) [arXiv :1109.3457], 

[17] E. Calabrese, R.A. Hlozek, N. Battaglia, et al., Phys. Rev. D , 87, 10 (2013) [arXiv: 1302.1841]. 

[18] M. Tegmark, M.R. Blanton, M.A. Strauss, etal.,ApJ, 606, 2 (2004) [astro-ph/0310725]. 

[19] C. Howlett, A. Ross, L. Samushia, et al., arXiv:1409.3238 (2014). 

[20] M. Levi, C. Bebek, T. Beers, et al, arXiv:1308.0847 (2013). 

[21] R. Laureijs, J. Amiaux, S. Arduini, et al., arXiv:1110.3193 (2011). 

[22] M. Bilicki, T.H. Jarrett, J.A. Peacock, et al., ApJS, 210, 9 (2014) [arXiv: 1311.5246]. 


6 



